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0\ ■ Abstract 

An analysis of the dilepton measurements of the reactions Pb(158 A-GeV) + Au, Pb 
^ ! by the CERES and NA50 collaborations points to a unique thermal source contributing 

in 
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to the invariant mass and transverse momentum spectra. 
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Qh! Introduction: Dileptons are penetrating probes which carry nearly undisturbed in- 

^ ■ formation about early, hot and dense matter stages in relativistic heavy-ion collisions. 
Some effort, however, is needed for disentangling the various sources contributing to the 
observed yields and for identifying the messengers from primordial states of strongly in- 
teracting matter. 

The dielectron spectra for the reaction Pb(158 A-GeV) + Au measured by the CERES 
collaboration ||1[ cannot be described by a superposition of e~^e~ decay channels of final 
hadrons, i.e. the hadronic cocktail. A significant additional source of dielectrons must 
be there. Since the data fll cover mainly the invariant mass range M < 1.5 GeV the 
downward extrapolation of the Drell-Yan process is not an appropriate explanation. Also 
correlated semileptonic decays of open charm mesons have been excluded . As a widely 
accepted explanation, a thermal source is found to account for the data (cf. |^ |] and 
further references therein, in particular with respect to in-medium effects and chiral sym- 
metry restoration). 

Very similar, the NA50 collaboration has found, for the reaction Pb(158 A-GeV) + Pb, 
that the superposition of Drell-Yan dimuons and open charm decays does not explain the 
data in the invariant mass range 1.5 GeV < M < 2.5 GeV 1^. Final state interactions 
or abnormal charm enhancement have been proposed as possible explanations. Here 



we try to explain the NA50 measurements by another idea p|, namely a thermal source 
|TU| , which was already found to account for the data in the intermediate invariant mass 



range in the reaction S (200 A-GeV) + W |TT|. We present a schematic model describing 



at the same time the CERES and NA50 data. 

The model: Since we include the corresponding detector acceptances a good starting 
point for Monte Carlo simulations is the differential dilepton spectrum 

dN r . dR 



P±idp±ip±2dp±2dyidy2d(l)id(j)2 J d^Q d^x 

where Q = P1+P2 is the pair four-momentum, pi 2 are the individual lepton four-momenta 
composed of transverse momenta p± 1^2 and rapidities 1/1^2 and azimuthal angles 0i 2- Here 
we extensively employ the quark - hadron duality |TT], and base the rate R on 

the lowest-order quark - antiquark (qq) annihilation rate (cf. [jl3|, |14|) 

where u{x) is the four- velocity of the medium depending on the space-time as also the 
temperature T{x) does. Note that, due to Lorentz invariance, u necessarily enters this 
expression. The above rate is in Boltzmann approximation, and a term related to the 
chemical potential is suppressed. As shown in the qq rate deviates from the hadronic 
one at M < 300 MeV, but in this range the Dalitz decays dominate anyhow; in addition, 
in this range the thermal yield is strongly suppressed by the CERES acceptance. In the 
kinematical regions we consider below, the lepton masses can be neglected. 

Performing the space-time and momentum integrations in eqs. (|l|, H) one gets 

dN 5a^ r^f 



PliP±2 f ' dtV{t)E, (3) 



dp± 1 dp± 2dyi dy2 dcpi d(j)2 7271^ Jt, 

{ f Mx cosh y cosh p(r,t) 1 sinhg r o f) 

TV- / M± cosh p{r,t) \ J ( Q± sinh p(r,t) \ r r) p, ^ ^ 

{ 4:71 f drr'^ for 3D, , , 

\t\drr for 2D, 

where V[t) acts on E, and 3D means spherical symmetric expansion, while 2D denotes 
the case of longitudinal boost-invariant and cylinder-symmetrical transverse expansion; 



the quantity ^ is defined as ^ = T ^ sinh p\jM'\_ cosh Y — M^, and p(r, t) is the radial or 



transverse expansion rapidity; Kq and Jq are Bessel functions ||13[ . The components of the 
lepton pair four-momentum Q = {M± cosh Y, M± sinh Y, Q±) are related to the individual 
lepton momenta via 

Ml = pii +J9I2 + 2j9±iP±2COsh(?/i - t/a), (6) 



Q± = P±l+P±2, (7) 

= Ml-Ql, (8) 

p±i cosh Hi + p^2 cosh ?/2 

It turns out that the shape of the invariant mass spectrum dN/{dM dY\\Y\<o.5dtdV{t)), 
which is determined only by the emissivity function E, does not depend on the flow 
rapidity p in the 2D case |T^, and in the 3D case for T = 120 • • • 220 MeV and p < 0.6 



there is also no effect of the flow. The analysis of transverse momentum spectra of various 
hadrons species point to an average transverse expansion velocity v± ~ 0.43 at kinetic 
freeze-out, while a combined analysis of hadron spectra including HBT data yields v± ~ 



0.55 [0. Therefore, p < 0.6 is the relevant range for the considered reactions. 

We note further that the invariant mass spectra dN/ {dM dY\\Y\<o.5 dt dV{t)) for the 
3D and 2D cases differ only marginally. Relying on these findings one can approximate 
the emissivity function E by that of a "static" source at midrapidity, as appropriate only 
for symmetric systems, 

f M^coshFl , , 

thus getting rid of the peculiarities of the flow pattern. We would like to emphasize the 
approximate character of eq. (p!0D because once cooling and dilution of the matter are 
included, they are necessarily accompanied by expansion and flow. One has therefore to 
check to which degree the flow affects the dilepton observables. 

In contrast to the invariant mass spectra, the transverse momentum or transverse 



mass spectra are sensitive to the flow pattern |T^, |17[ |18[, in general. A value of 
p = 0.4, for example, causes already a sizeable change of the shape of the spectra 
dN/{dQ±dY\\y\^Q^^dt dV{t)) compared to p = 0, in particular in the large-Qx region. 
The differences between the 2D and 3D cases are not larger than a factor of 2 and, in a 
restricted Q± interval, can be absorbed in a renormalization. The most striking difference 
of the 2D and 3D scenarios is seen in the rapidity spectrum: for 2D it is fiat, while in 
the 3D case it is localized at midrapidity (values of p < 0.6 also do not change the lat- 
ter rapidity distribution). Below we shall discuss which space is left to extract from the 
dilepton data in restricted acceptance regions hints for the flow pattern when the other 
dilepton sources are also taken into account. 

Comparison with data: In line with the above arguments we base our rate calcu- 
lations on eqs. (j^, [1^) and use the parameterizations 



T = {Ti-T^)exp\-l\+T^, (11) 



t 



V = A^expj^j. (12) 
with Ti = 210 MeV, Too = HO MeV, h = 5 fm/c, = 8 fm/c, = |g with A, B as 



mass numbers of the colliding systems and Uq = 0.17 fm ^. We stop the time evolution 
at Tf = 130 MeV. 

In fig. 1 we show the comparison with the preliminary CERES data applying the 
appropriate acceptance |jll|. One observes a satisfactory overall agreement of the sum 
of the hadronic cocktail and the thermal contribution with the data. It is the thermal 
contribution which fills the hole of the cocktail around M = 0.5 GeV in the invariant mass 
distribution in fig. la. In the mass bin M = 0.25 ■ ■ ■ 0.68 GeV the thermal yield is seen 
(fig. lb) to dominate at small values of the transverse momentum Q±. In this region of 
Q± transverse fiow effects are not important. The large-Q_L spectrum is dominated by the 
cocktail. For higher-mass bins the thermal yield in the region of the first two data points 
is nearly as strong as the cocktail and rapidly falls then at larger values of Q± below the 
cocktail. Therefore, the fiow effects turn out to be of minor importance for the present 
analysis, since within our framework the transverse flow shows up at larger values of Q±. 

The question now is whether the same thermal source model accounts also for the 
NA50 data In the mass range M > 1 GeV, the Drell-Yan dileptons and dileptons 
from correlated semileptonic decays of open charm mesons must be included. To get the 



corresponding yields for Pb + Pb collisions from PYTHIA the overlap function T^a = 
31 mb~^ is used. We have carefully checked that our PYTHIA calculations with K factors 



i^DY = 1-23 and -ft'charm = 4, adjustcd to Drell-Yan data and identifled open charm 
data (cf. [0 for a data compilation), and intrinsic transverse parton momentum {k±) = 0.9 
GeV coincide with results of the NA50 collaboration. In particular, we reproduce with 
Xd.o.f. — 0-24 the anticipated NA50 result that the data of the most central collisions 
(except the J/ip region) are described by the Drell-Yan yield + 2.8 x the yield from 
correlated semileptonic decays of open charm. In this way we get some confldence in 
our acceptance routine which essentially consists of geometrical cuts and a suitable 
minimum single-muon energy of (9(11.5 GeV) [0. 

The resulting invariant mass spectra, including the thermal source contribution, are 
displayed in flg. 2a. The thermal source, with strength adjusted by the above comparison 
with CERES data, is needed to achieve the overall agreement with data. For M < 2 
GeV the thermal contribution dominates over the Dell-Yan and charm contributions. 
The value of Xdof — 1-38 quantifles that the very details of the data are not perfectly 
described. This may be attributed to the too schematic source model eq. (|TOD and our 
approximate description of the more involved NA50 acceptance. Nevertheless, the trans- 
verse momentum spectrum for the mass bin M = 1.5 ■ ■ ■ 2.5 GeV is nicely reproduced 
by the sum of Drell-Yan, open charm and thermal contributions with Xdof =2. The 
thermal yield is strongest at not too large values of Q± where transverse flow effects can 
be neglected. Therefore, it seems that from present dilepton measurements the transverse 
flow can hardly be inferred. However, a reduction of the above uncomfortably large value 



of the phenomenological parameter {k±) could be partially compensated by transverse 
flow. 

In a previous version of this note ||^ we have confronted our model with the flndings of 
the NA50 collaboration according to the recipe "reconstructed data" = the Drell-Yan yield 
+ 3 X the dilepton yield from decays of D mesons and found a much better agreement. 
This fact points to the need to employ the full acceptance matrix in attempting a decision 
on the best data interpretation. 

Summary and discussion: In summary we have shown that a very simplifled, 
schematic model for thermal dilepton emission, with parameters adjusted to the CERES 
data, also accounts for the measurements of the NA50 collaboration. Our study points to a 
common thermal source seen in different phase space regions in the dielectron and dimuon 
channels. This unifying interpretation of different measurements has to be contrasted 
with other proposals of explaining the dimuon excess in the invariant mass region 1.5 ■ ■ • 
2.5 GeV either by flnal state hadronic interactions or an abnormally large open charm 
production. The latter one should be checked experimentally by a direct measurement of 
D mesons as envisaged in the proposal |]^, thus providing a flrm understanding of various 
dilepton sources. 

Due to the convolution of the local matter emissivity and the space-time history of the 
whole matter and the general dependence on the flow pattern, it is difficult to decide which 
type of matter (deconflned or hadron matter) emits really the dileptons. Our model is not 
aimed at answering this question. Instead, with respect to chiral symmetry restoration, 
we apply the quark - hadron duality as a convenient way to roughly describe the dilepton 
emissivity of matter by a qq rate, being aware that higher-order QCD processes change 



this rate |]T8|, ^ (this might partially be included in a changed normalization A^). In 
further investigations a more microscopically founded rate together with a more detailed 
space-time evolution must be attempted and chemical potentials controlling the baryon 



and pion densities must be included. This has been accomplished in to a large extent 
with similar conclusions as ours. 
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Figure 1: The preliminary CERES data [0] and the hadronic cocktail [|l| (dashed lines) 
and the thermal yield (full curves). The sum of the cocktail and the thermal yield is 
shown by the gray curves. Left panel (a): the invariant mass spectrum (Xdof — l-^), 
right panel (b): the transverse momentum spectrum for the mass bin 0.25 • • ■ 0.68 GeV 
(Xio.f. = -74). 
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Figure 2: The preliminary NA50 data [^] in comparison with the thermal yield (full 
curves), the Drell-Yan contribution (dashed curves) and the contribution of open charm 
decays (dash-dotted curves). The J/tp and ip' curves (thin lines) are taken from [^]. The 
sum of these contributions is displayed by the gray curves. Left panel (a): the continuum 
invariant mass spectrum, right panel (b): the transverse momentum spectrum for the 
mass bin 1.5 ■ ■ • 2.5 GeV. 



